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Exploiting isochrony in self-triggered control 

Adolfo Anta Paulo Tabuada 



Abstract 

Event-triggered control and self-triggered control have been recently proposed as new implementa- 
tion paradigms that reduce resource usage for control systems. In self-triggered control, the controller 
is augmented with the computation of the next time instant at which the feedback control law is to be 
recomputed. Since these execution instants are obtained as a function of the plant state, we effectively 
close the loop only when it is required to maintain the desired performance, thereby greatly reducing 
the resources required for control. In this paper we present a new technique for the computation of the 
execution instants by exploiting the concept of isochronous manifolds, also introduced in this paper. 
While our previous results showed how homogeneity can be used to compute the execution instants 
along some directions in the state space, the concept of isochrony allows us to compute the executions 
instants along every direction in the state space. Moreover, we also show in this paper how to homogenize 
smooth control systems thus making our results applicable to any smooth control system. The benefits 
of the proposed approach with respect to existing techniques are analyzed in two examples. 

I. Introduction 

Recent technological and economical trends require control systems not to be implemented 
on dedicated platforms, but rather on shared devices where different subsystems compete for the 
available resources (e.g. computation time, communication bandwidth, shared actuators). In this 
shared setup, traditional implementation paradigms are no longer adequate since they require 
the same amount of resources independently of the workload, performance, or the effect of 
exogenous disturbances. For instance, in a digital implementation, control laws are executed 
at the same rate regardless of the current state of the system. While this approach simplifies 
the analysis and implementation of the control system, the choice of the period for the control 
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loop is independent of the current state and performance of the system. To overcome the short- 
comings of periodic implementations, several researchers advocated the use of event-triggered 
control IAiz99j . IAB021 . HTab07ll . HHSvdBOSl . [HXTOTl . Under this paradigm, control tasks are 
executed whenever a triggering condition is satisfied (e.g. whenever the state deviates from a 
nominal trajectory by a certain amount). The use of event-triggered control, however, entails 
one handicap: dedicated hardware is required to continuously check the triggering condition. 
Even if such hardware can be deployed, it implies a significant cost. To surmount this drawback, 
Velasco et al [VFM03| suggested the self-triggered control concept, by which the next execution 
time is given by a function of the last measurement of the state. Self-triggered control can 
be regarded as a software-based emulation of event-triggered control. Self-triggered control for 
linear systems was studied in HWL09H , HMJAT10L It has also been recently applied in the context 
of sensor/actuator networks HMJT08L HTFJDB10L HCMV+10H . To the best of our knowledge, 
self-triggered formulas for nonlinear systems were first developed in HAT 101 . 

In this paper we extend our previous work on self-triggered control in several directions. 
In HAT 1 PL the geometric properties of homogeneous vector fields were exploited to derive 
scaling laws describing the inter-execution times, that is, the times elapsed between consecutive 
executions of the control law. These results explained how such times evolve along homogeneous 
rays. In this paper we develop formulas that explain how the inter-execution times evolve across 
homogeneous rays. The contributions of this paper are as follows: 

• We introduce the concept of isochronous manifolds, show their existence, and how it can 
be used to describe how the inter-execution times change across homogeneous rays. 

• We combine our previous results in [AT10] with the concept of isochronous manifolds to 
obtain new self-triggering formulas that explain how the inter-execution times evolve along 
rays and across rays. 

• We generalize the results in HAT 1 Oil to non-homogeneous control systems and non-homogeneous 
triggering conditions. 

We also analyze the accuracy of the proposed self-triggering formulas by: 

• Providing a bound for the mismatch between the inter-execution times generated by the 
event-triggering technique and by the self-triggering technique. 

• Analyzing the inherent trade-off between complexity and accuracy of the self-triggering 
formula, and deriving an iterative algorithm that trades complexity for precision. 

• Discussing the role played by the different design parameters in the accuracy of the self- 
triggering formulas. 

This paper does not focus on the self-triggered emulation of a particular event-triggered 
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implementation, but rather develops a self-triggering technique that can be used with any smooth 
event-triggering condition defined for any closed-loop system described by a smooth vector field. 
We illustrate the results on two examples, showing that the proposed technique outperforms 
previous solutions and produces self-triggered times that hardly differ from the event-triggered 
times. A preliminary version of these results appeared in [AT09J. 

II. Notation 

We use the notation \x\ to denote the Euclidean norm of an element x E W 1 . A function is 
said to be of class C°° or smooth if it can be differentiated infinitely many times. All the objects 
in this paper are considered to be smooth unless otherwise stated. 

Given vector fields X and Y in an n-dimensional manifold M, we let [X, Y] denote their Lie 
product which, in local coordinates x = (xi,x 2 , x n ), we take as ^(x)X(x) — ^(x)Y(x). 
We use Cxh(x) to denote the Lie derivative of a map h : M — > R at a point x along the flow of 
the vector field X : M — > TM which, in local coordinates, we take as ||X(x). Likewise, C k x h 
represents the A;th Lie derivative, defined by C° x h = h and C x h = Cx(Cxh). Let : M — >• N 
be a map, and let X and Y be vector fields in manifolds M and N, respectively. We denote 
the differential of a map <p by T<f>, which, in local coordinates, we take as 1^. We call X and 
Y 0-related if the following holds: 

T<f>- X = Y o0. (Il l) 

In this paper we consider control systems of the form: 

x = f(x, u), x(t) g R n , u(t) eR m ,te Rj. (II.2) 



We denote by i G M" the state of the control system, by x a solution of ( |IL2[ ), and by u the 



input trajectory. Finally, we use e A to represent the exponential of a matrix A £ ]R nxn , Ay to 



denote the element ij of the matrix A and I to represent the identity matrix. 



III. Self-triggering stabilization of nonlinear systems 
We start by analyzing the behavior of a control system under a sample-and-hold implementa- 



tion. Consider the system (II. 2), for which a feedback control law of the form u = k(x) has been 
designed. The implementation of such control law on an embedded processor is typically done 
by sampling the current state of the system x(£) at time instants t — ti, t 2 > • • •> computing the new 
value for the input u(tj) = k(x(ti)), and updating the actuator with u(U) at time instants U. We 
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shall use the term execution to denote this process of sampling the state, computing the control 
law and updating the input. Between actuator updates the input u is held constant, according to: 

te [ti,t i+1 [=*u(t)=u(ti). (HIT) 

Traditionally, the controller is executed periodically, that is t i+1 —ti — T, for any i G N, where 
T > is the period. Instead, we would like to identify a sequence of execution time instants 
that guarantees stability and desired performance while reducing the number of executions. To 
proceed, we define the measurement error e as the difference between the last measurement of 
the state and the current value of the state: 

e(t) = x(ti) - x(t) for t E [U, t i+1 [. (HI.2) 

With this definition, the closed loop x = /(x, k(x(ti))) becomes: 

x = /(x,Jfe(x + e)). (HI.3) 



Whenever a control law u = k(x) is designed to render the system (II. 2) asymptotically stable, 
there exist a Lyapunov function V satisfying: 

V\ t=0 = ^f(x, k(x)) <0 VxG R n \{0}. (III.4) 

Moreover, control laws are usually designed not only to guarantee stability but also to achieve a 
certain performance. In this paper the performance requirement will be given by a desired rate 
of decay h : IR n — > E of the Lyapunov function. Hence V satisfies: 

dV 

— f(x, k{x)) < h{x) < 0, Wxe K n \{0}. (III.5) 



However, under a digital implementation the dynamics of the system are given by (III.3 ). Let the 
desired rate of decay for the implemented closed loop be h(x) > h(x). If the error e satisfies: 

dV 

— fix, k(x + e)) < h(x) < (III.6) 
ox 

the rate of decay h can be guaranteed for the s ample- and-hold implementation. Therefore, the 
sequence of time instants U has to be chosen so that inequality ( |IIL6 ) is satisfied at any time. 
This can be achieved by closing the loop whenever: 

dV 

T(x, e) := -7— fix, kix + e)) - hix) = 0. (III.7) 
ox 



Equation (III.7) implicitly defines the time instants r at which the state x needs to be measured, 
the control law u = k(x) computed, and the actuator updated with the control input u: 

r(x(U),e(U)) := min{t > U : r(x(t,x(t i ),e(t i )),e(t,x(t l ),e(t i ))) = 0}. 
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Whenever there does not exist t such that r(x(t, x(ij), e(tj)), e(t, x(^), e(^))) = 0, we will say 



that r(x(tj), e(U)) = oo. Indeed, if (III.7) is never satisfied the control law does not need to be 



updated anymore. Notice that our analysis is not constrained by the particular structure of ( |IIL7| ) 
but extends to other typical (smooth) triggering conditions available in the literature, as it will 
be shown in the remainder of the paper. 



Equation (III.7) represents an event-triggering condition that renders the system asymptot- 
ically stable and guarantees the desired rate of decay h. Upon the execution of the control 
law at t = U, the state is measured and the error becomes 0, since x(t) = x(tj) implies 
e(ti) = x(tj) — x(ij) = 0. An event-triggered implementation based on this equality would re- 



quire testing (III.6) frequently. Unless this testing process is implemented in hardware, one 



might run the risk of consuming the processor time testing (III.6). To overcome this drawback, we 
consider self-triggering strategies, where the current state measurement x(U) is used to determine 
the next execution time for the control law. Self-triggered implementations represent a software- 
based emulation of event-triggered implementations. 

For clarity of exposition, at this point we define a new state z E M 2n that includes the state 
of the dynamical system x and the measurement error e: 

x 
e 

Likewise, we use z to denote the trajectories of the extended state z. 

The self-triggering strategy decides the next inter-execution time according to a function 
T l . j^2n _j, j£+ fast lower bounds the event-triggered times, that is, the inter-execution times 
of the event-triggering policy: 



(III.8) 



where z(U) = (x(ti),e(ti)) corresponds to the last measurement of the current state x(U) and to 
the error at the sampling time e(ij). Since 

r(z(t,z(t<))) < Wte[U,r[ (III.9) 

closing the loop at r, or before, enforces ( |IIL6[ ) for all t > 0. 

Therefore, the self-triggered times should be as large as possible (in order to reduce the 
number of executions of the control law) but no larger than the event-triggered times (in order 
to guarantee the desired performance). Since the state trajectories are in general not known for 
nonlinear systems, it is not possible to compute the inter-execution times r(z(tj)) in closed form. 
In the next section we develop scaling laws for the inter-execution times of smooth nonlinear 



systems implicitly defined by the triggering condition ( |III.7| ). 
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IV. Scaling laws for the inter-execution times 

We start by reviewing the scaling laws previously derived in HAT 1011 for the special class of 
homogeneous systems. 



A. Homogeneous systems 

We first define homogeneity for functions. 

Definition IV. 1. A function f : M. n — > R n is called homogeneous of degree £ G R if there exist 
Tn > 0, i — 1 . . . n such that: 



VA > f i (X n x 1 , Y"x n ) = A^/^i, ...,x n ) 
where £ > — min, rj. 



(IV.l) 



In a geometric setting, a function / is said to be homogeneous of degree £ if the following 
holds: 

C D f = (£ + l)f (IV.2) 

where D is known as a dilation vector field, i.e., a vector field such that x = —D(x) is globally 
asymptotically stable. When the choice 

d 



»=i 



dxi 



(IV.3) 



is made, the geometric definition in ( |IV.2[ ) results in Definition |IV. 1 1 (see HKaw95B , [|Arn92| for 
details). Likewise, we define the notion of homogeneity for vector fields. 



Definition IV.2. hKaw95\l Let D : M — >■ TM be a dilation vector field. A vector field 
X : M — > TM is called homogeneous of degree £ G R with respect to the vector field D if 
it satisfies the following relation: 

[D,X]=£X. (IV.4) 
To derive scaling laws for the self-triggered times, we consider the standard dilation vector 



field (IV.3). The trajectory of the vector field D starting at the initial condition x G 



corresponds to the standard open ray {Xxo : A G IR + }. Using the inherent symmetries of 
homogeneous vector fields, a scaling law was developed in HAT 1 OH under the following triggering 
condition: 

" -a 2 I 



r(z(r, zq)) = z 1 (t, z )Mz(t, zq) = 0, 



M 







(IV.5) 
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where a G IR + represents a design parameter trading control performance for number of execu- 
tions. 

Theorem IV.3. [AT10] Let z = Z(z) be a dynamical system homogeneous of degree £ G K 



with respect to the standard dilation vector field. The inter-execution time t 



)2;i 



— > 



implicitly defined by the execution rule ( |IV.5[ ) scales according to: 

r(Az ) = A-V(zo), A > 
where Zq G M 2n represents any point in the extended state space. 



U {00} 



(IV.6) 



In the remainder of this section we extend these results for non-homogeneous systems and 
non-homogeneous triggering conditions. 

B. Homogenization of dynamical systems 

We start describing how smooth nonlinear systems can be rendered homogeneous by intro- 
ducing an auxiliary state variable. The procedure is well known for polynomial systems (see for 
example [Bai80]). For instance, the following system: 



Xl 

x 2 



XiX 2 + x 2 



(IV.7) 



can be rendered homogeneous of degree 1 with respect to the standard dilation by adding an 
auxiliary variable weR satisfying w = 0: 



Xl 

x 2 
w 



x x x 2 + X 2 W 
XiW 

0. 



(IV.8) 



Trajectories for the original system ( IV.7[ ) can be recovered from the trajectories of the auxiliary 



system (IV.8) provided that w(0) = 1. Formally speaking, we say that the vector field 

X -- 



XiX 2 + x 2 



(IV.9) 



is (^-related with the vector field 



S 



XiX 2 + x 2 w 
XiW 





(IV. 10) 
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with (p(x) = (x, 1), since Tip ■ X = S o (p. This idea was exploited in our previous work HAT 1 OB 
to derive scaling laws for polynomial systems using Theorem IV.3 We generalize this homog- 
enization procedure to smooth nonlinear systems in the following lemma[j] 



Lemma IV.4. For any vector field 



i=l 



I — 

dxi 



there exists a vector field S ip-related to X, with ip(x) = (x, 1), homogeneous of degree £ G 
with respect to the standard dilation vector field, and given by: 

n d 
S(x,w) = w^ +1 fj{w~ l x) 



i=i 



dxi 



(IV.ll) 



Proof: We first prove that S is homogeneous of degree £. The standard dilation vector field 
for S is given by: 

^— ' OXa ow 



i=l 



We compute the Lie bracket between D and S to verify condition ( |IV.4 ): 

[D,S] 



dx dx 



£ x (w^f(w-*x)) £(w^f(w^x)) 





x 




w 







s 



d(w x) 

since -Sj — 





d(w~ 1 x) 





5 = £S 



s 



(IV. 13) 



<9x 



satisfied: 



; w. Moreover, the vector fields X are S are (^-related since <|II. 1 p is 

^ ^ = S(x, 1) = S o cp. 



T<p-X 



(IV. 14) 



It is easy to see that the homogenization procedure for polynomial systems represents a 



particular case of Lemma IV.4 Notice that the auxiliary system is not defined for w — 0. 
Nonetheless, it will be shown in section IV] that the 2n + 1 dimensional model is only used 



Although we believe that such result might have appeared in the literature, we were not able to find it. For completeness, 
we include a proof. 
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for points (z, w) with w ^ 0. By embedding the original 2n-dimensional dynamical system 
into a 2n + 1 dimensional homogeneous system we can derive similar scaling laws for the 
inter-execution times, as explained in the next subsection. 

C. Homogenization of triggering conditions 

We also extend the class of triggering conditions to be considered. In our previous work HAT 1 Oil 
we introduced the simple event- triggering condition ( |IV.5| ). However, in general one might require 



a different triggering condition. For instance, another possible condition is V = 0, in order to 
guarantee V < for all t G IR + . 



Theorem IV.3 leverages the fact that the triggering condition ( |IV.5[ ) is quadratic in z. Consider 
now a triggering condition of the form T(z(t, z)) = 0, where T is a homogeneous function of 
degree and the trajectory z is governed by a vector field homogeneous of degree £. Then, the 
scaling laws for the inter-execution times are the same regardless of the degree of homogeneity 
of the function T. Indeed, note that the triggering condition T satisfies: 

r(z(t, Xz )) = T(Xz(Xk, zo)) = A*T(z(A*t, z )) (TV. 15) 

where the first equality is a property of homogeneous flows (see [Kaw95] for details). Since 
X s 7^ 0, whenever 

X^T{z{Xk,z )) =0 (IV. 16) 

holds, the triggering condition 

r(z(t,Az )) =0 (TV.17) 
is fulfilled, and vice versa. In other words, the time instants XH at which (IIV. 1 61) is satisfied 



coincide with the time instants t at which ( |IV17| ) is satisfied. Hence, the same scaling laws hold 



whenever the triggering condition is a homogeneous function. 

Moreover, a similar procedure as the one used for vector fields can be exploited to render 
a more general class of triggering conditions homogeneous. Starting with a nonhomogeneous 
triggering condition T(z) = we can construct a homogeneous triggering condition of degree 
•d by means of an auxiliary variable w: 

w^Tiw^z) = (IV.18) 

that is equivalent to the original triggering condition T(z) = for w = 1. By working in a 2n + l 
dimensional space we are able to construct symmetries for the inter-execution times. These two 
homogenization procedures for the triggering condition and the vector field are summarized in 
the following theorem. 
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Theorem IV.5. Let r : M 2n — > IR + U {00} be the inter-execution times for the system z = Z(z) 
implicitly defined by T(z(t(z), z)) = for a point z G R 2n , and let f : R 2n+l — > R + U {00} be 
the inter-execution times for the system: 



z 
w 



S(z, w) 



w« +1 Z(w- 1 z) 




(IV. 19) 



implicitly defined by T(z,w) := T(w 1 z(f(z, w), z, w)) = 0, for a point (z,w) G IR 2n+1 . Then, 
r and f are related according to: 



-( z ) = \tf(\ z ,\) VA>0. 



(IV.20) 



Proof: We first note that Z : M -> TM is (^-related to S : N TN with = (z, 1). 
Thus the following diagram commutes: 

r Z 



(IV.21) 



since the following two conditions are satisfied: 

1) Tip ■ Z — S o (p. 

2) ro<p = r. 

That is, in order to study the evolution of the map T(z) under the vector field Z we can 
alternatively study the evolution of the map T(w~ 1 z) under the dynamics of S. Let s be the flow 
for the vector field S. Since Z and S are (^-related, the corresponding flows satisfy [Lee03J: 




ip o z(t, z) = s(t, (f(z)), z G M 



and from commutativity of Diagram IV.21 we can conclude: 



r o z(t, z) = r o tp o z(t, z) = r o s(t, <p(z)). 



(IV.22) 



(IV.23) 



Since the evolution of both triggering conditions T o z and f o s o ip is identical, they generate 
the same inter-execution time for any point z G M: 

t{z)=t(z,1). (IV.24) 

Finally, since the vector field S is homogeneous of degree £, times scale according to Equa- 



tion (IV.6): 



-(z) — f(z, 1) = \tf(\z, A) VAX). 



(IV.25) 
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Notice that this result does not rely on the fact that e = — x, but rather on the fact that the 
augmented vector field S is homogeneous. In the following example we show how the previous 
theorem can be applied to study the evolution of the trajectories of a nonlinear system. 

Example IV.6. Consider a simple unforced pendulum described by: 

±i = x 2 

x 2 = -7 sin (xi) - — x 2 (IV.26) 
I m 

where x\ corresponds to the angle of the pendulum (with respect to a vertical line), x 2 its angular 
velocity, m is the mass of the pendulum, I the length of the rod, k is the coefficient of friction 
and g the gravitational constant. We are interested in analyzing the time at which the angle of 
the pendulum reaches 71/6, starting with initial condition (xi ,x 2o ). Hence, we are looking for 
the time r at which the following equality is satisfied: 

T(x 1 (t,x 1o ,x 2o )) := Xi(r,Xi ,x 2o ) - tt/6 = 0. (IV.27) 

We homogenize the system as described before (for instance, with degree of homogeneity £ = 1): 

Xi = X 2 W 

9 / — 1 \ 2 r\ r» 

x 2 = — -sinlw xijw x 2 w (IV.28) 

w = 0. 



Likewise, equality (IV.27 ) becomes: 



~ ~ 71" 

r(xi(r,Xi ,a; 2o ,Wo)) := xi(r, x lo , x 2o , w ) - -w(r, x lo , x 2o , w ) = 0. (IV.29) 

o 

The time f implicitly defined by r(xi(f, X\ , x 2o )) = and the time r implicitly defined by 
r(xi(r, X\ , x 2o )) = are related according to: 

t(x 1q ,x 2o ) = Xf(Xx lo , Xx 2o , A) A > 0. (IV.30) 

As shown in this example, the applicability of the previous framework goes beyond the self- 
triggered control problem and can be used to describe other temporal aspects of the trajectories 



of nonlinear systems. In the next section we describe how equation (IV.30) can be used to 



determine the time at which equality (IV.27 ) is satisfied 
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V. Isochronous manifolds 
In this section we exploit the results of the previous section to derive a self-triggering strategy 



that emulates an event-triggered implementation. In Section IV scaling laws were derived for 
the inter-execution times of dynamical systems. Such scaling laws describe the evolution of the 
times implicitly defined by a triggering condition. Once the inter-execution time for a point is 
known, times for any point lying on the same ray can be derived by means of the scaling laws. 
Hence, in order to compute the inter-execution times for the whole extended state space it is 
necessary to obtain the inter-execution times for at least one point in each ray. For instance, it 
would be enough to find a lower bound r* for the times on a sphere (since every ray intersects 
a sphere), and then extend the times along rays. This was the approach used in our previous 
work HAT lOH . Different methods are available to compute a lower bound r* for the inter-execution 
times on a sphere (see for instance HTab071l ). Following this procedure, inter-execution times for 



the self-triggering strategy can be computed using (IV.6): 



t(z(U)) = A~V\ with A = (V.l) 



where we have identified the point Xzq in (IV.6) with the last measurement of the extended state 



z(ij) = (x(tj), 0), the left hand side of (IV.6) with the time for z(t$), z with a point in a sphere 
of radius r, and t(zq) with r*, the lower bound for such sphere. Unfortunately, such procedure 
introduces two sources of conservativeness: 

• The existing methods to compute a stabilizing period for a given set (also known in the 
literature as maximum allowable transfer interval, MATI) are rather conservative HNTC 09I, 
[Tab07], since in general they are based on the Lipschitz constant of the vector field. 

• The stabilizing period represents a lower bound for the inter-execution times of all the points 
in the sphere. Since times might vary drastically along the sphere (or along any given set), 
this lower bound will be overly conservative for some points in the sphere. In other words, 
the evolution of the times across rays is neglected in this procedure. 

In the next subsections we overcome these two drawbacks through the use of isochronous 
manifolds. 



A. Existence of isochronous manifolds 

As mentioned before, the sphere was selected in HAT 1011 as the reference manifold (i.e., the set 
where a lower bound for the inter-execution times is computed) for simplicity. However, there is 
no reason to expect that such times are constant for the points in the sphere. Ideally, we would 
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like to identify manifolds where the inter-execution times remain constant. Such manifolds will 
be termed as isochronous manifolds. If these manifolds can be computed in an exact manner, 
the inter-execution times generated by the self-triggering strategy will coincide with the event- 
triggered times. 

It is important to notice that we are interested in isochronous manifolds that allow us to derive 



self-triggering conditions for the whole operating region through the scaling laws of Section IV 
In other words, in order to be of any use, an isochronous manifold should intersect every ray at 
least at one point. We prove the existence of this class of isochronous manifolds for homogeneous 
systems in the following proposition. 

Proposition V.l. Consider the dynamical system z = Z(z) homogeneous of degree £ > 0, 
defined on a 2n- dimensional manifold. Let r : IR 2n — > IR + U {00} be implicitly defined by 
T(z(t(zo), Zo)) = 0, where T is a homogeneous function. If the following two assumptions are 
satisfied for any Zo in a sphere of radius r (for some r G M. + ): 

r(z(0,z )) t^O, (V.2) 

3t G IR + : T(z(t, zq)) = (V.3) 

then, for any given time t t . G M + the set Q tt defined by: 

ft u = {z, G R 2n : t(z,) = U} (V.4) 

is a manifold of dimension 2n — 1. 



Proof: As explained in Section IV-C the inter-execution time r defined by the triggering 
condition scales according to: 

r(Az ) = X~ $ t{z ), VA > 0, f > 0. (V.5) 



According to (V.5), inter-execution times vary from to 00 as A varies from 00 to 0. Hence, 



for any given time £* G M + , there exists a point z* in each ray such that t(z*) = t*. Moreover, 



equation (V.5 ) implies that there does not exist two different points lying on the same ray with the 
same inter-execution time. The union of all those points z* constitutes an isochronous manifold 
for the homogeneous system. 

We briefly sketch now the proof that fl, t * is indeed a manifold. We use d to denote the exterior 
derivative of a function. We first observe that O t „ can be equivalently defined as the zero level 
set of the map T(z(t*, z)) for the chosen t*. Since V is homogeneous we have: 

T(Xz) = X^T(z), VA > (V.6) 
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where $ is the degree of homogeneity of T. Then, ( |V.2[ ) implies that the triggering condition 
never vanishes for any z^O: 



r 



\z\rz 
r \z\ 



r 



r2 

Ul 



7^0 



(V.7) 



as the point f4 belongs to a sphere of radius r. Since T is a homogeneous function it satisfies: 



c D v = (c + i)r 



(V.8) 



with D being the standard dilation vector field, as defined in ( |rV.3[ ). Since the Lie derivative can 
also be expressed as: 

dr(£>) = C D Y. (V.9) 

and D only vanishes at the origin, we conclude from ( |V.8[ ) and ( V.9| ) that dr only vanishes 
when T = 0, i.e., at z = 0; therefore, T is a submersion for every z ^ 0. Since z(t*,2) is a 
diffeomorphism for every t* E M, the corresponding tangent map is an isomorphism for every 
z E R n . Hence, we conclude that T(z(U,z)) is a smooth submersion from which follows that 
its zero level set is a smooth manifold of codimension 1 HLee031l . ■ 
Let us comment now on the two assumptions < |V.2| ) and ( |V.3[ ). Assumption < |V.2[ ) will be 
satisfied by the triggering conditions that are usually of interest. Indeed, should the triggering 
condition r(z(0, Zq)) be for some z , the inter-execution time for such z will be 0, which 



is undesirable from a practical point of view. For instance, for the triggering condition ( |IV.5| ) 
mentioned in Section [TV] we have e(0, Zq) 



r(z(o,zo)) 



0, hence assumption ( |V.2| ) holds: 

z(0, ^o) T Mz(0, zq) = e(0, ^ ) T e(0, z ) - ^(0, ^ ) T x(0, zq) 

-(T 2 X (0, ^ ) T X(0, Z ) < V2 7^ 0. 



(V.10) 



Likewise, it can be easily checked that the triggering condition in ( |III.7[ ) satisfies this condition. 

Assumption ( |V.3[ ) guarantees that there does not exist states where the inter-execution times 
are oo. In general, there might exist stable manifolds where no updates in the control input are 
needed to steer the system to the equilibrium point. In this case, the isochronous manifold would 
not intersect every homogeneous ray, and therefore it would not be useful to derive self-trigger 



formulas. For those states lying in stable manifolds, the scaling law (V.5) implies that times will 
be oo for all points lying on the rays intersecting the stable manifold. And vice versa, since 
all those points have oo as the next inter-execution time, any stable manifold for homogeneous 



systems needs to be composed of homogeneous rays. Thus, Assumption ( |V.3[ ) guarantees that 
times on each ray assume values in ]0, oo[. 



April 27, 2011 



DRAFT 



15 



The isochronous manifold can be used instead of the sphere in order to compute the inter- 
execution times for the whole operating region, as there exists a unique intersection point between 
each ray and the manifold. In the next section we provide techniques for the computation of 
such manifolds. 



B. Approximating isochronous manifolds 

The isochronous manifold Q tit implicitly defined by equation ( V.4[ ) represents the set of all 



extended states z* that satisfy the triggering condition at time £*. Alternatively, the isochronous 
manifold can be defined as: 

fit, = {z* eR n |r(z(t*,^)) = o}. (v.ii) 

The explicit computation of fit, would require the knowledge of the flow, which is in general 
unknown. Herein we develop a technique to approximate the isochronous manifold. Since we are 
only interested in the evolution of the triggering condition T(z(t, z )) = 0, it would be sufficient 
to construct a differential equation for T. For ease of exposition and without loss of generality, 
at this point we assume that T(z(0, z )) < for any z in the operating region (which is the case 
for the triggering conditions ( |IV.5[ ) and V = 0). In [Tab07] a first order nonlinear differential 
inequality of the form j^T < v(T) was derived, whose solution just depended on the norm of 
the initial extended state z . Instead, we derive a p-th order differential equation, whose solution 
collects more information from the extended state. Moreover, in order to obtain a closed-form 
expression for its solution, we focus on linear differential equations describing the evolution of 
T. Clearly, T will be described by a linear differential equation only under special circumstances, 
see e.g. HLM 86I. However, the evolution of T can be bounded by a linear differential equation. 
This amounts to constructing a set of coefficients Xi G IR satisfying: 

p-i 

c p z r(z) <^Xi4r^ (v.i2) 

with Z representing the closed-loop vector field and for any extended state z in an invariant 
region of interest M C ]R 2n , as stated in the next result. 

Lemma V.2. Consider a vector field Z : M — >■ TM and a map T : M — > R. For every set of 
coefficients Xo, Xi, ■ ■ ■ X P -i € R satisfying: 

p-i 

C p z T(z) < Y,XiC l z T{z) Vz E M, (V.13) 

i=0 
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the following inequality holds for all Zq G M and for all t G for which the solution z of Z 
is defined: 

r(z(t,zb)) < Yi(t,y ) 



where yi is the first component of the solution of the linear differential equation: 

Yi = Yi+i, z = l,...,p-l 

with initial condition: 

y = (T(z ),C z T(z ),...,C p z 1 T(z )) T . 



Yi+i, 
p-i 



i=0 



Proof: The differential inequality ( |V.13| ) can be rewritten as: 

U = bi+i z = l,...,p-l 

p-i 

bp < X ihi 

i=0 



(V.14) 

(V.15) 
(V.16) 

(V.17) 

(V.18) 
(V.19) 



with initial condition b = (T(zq), CzT(z ), . . . ,£ p z 1 r(2; )) T . By the Comparison Lemma (see 
for instance [Kha02]) we have b p (t, b ) < y p (t, y ) since b = y and b p (t, b ) = C p z T(z(t, z )) < 



SjXi^r(z(t,zo)) = y p (t,y ), in virtue of (V.13). We apply the Comparison Lemma succes- 
sively to conclude that each y« bounds each bf. 

bi(t,b ) <yi(t,yo), i = l,...,p, t>0. (V.20) 
Hence, yi represents a bound for the evolution of T: 

r(z(t, Zq)) = bi(t, b ) < Yi(t, yo) 



For convenience, at this point we define the following map: 



{T{z\CzT{z\...,a> z - l T{z)Y. 



(V.21) 



Using Lemma V.2 we can see that the minimum time satisfying yi(t,y ) = lower bounds the 
minimum time satisfying r(z(t,z )) = 0, since r(z(0, Zo)) < 0. Thus, the inter-execution time 
defined by yi(r^, n P (z)) = represent a lower bound for the inter-execution time r defined 
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by the triggering condition r(z(r, z )) = 0. We can interpret this lemma as a way to construct 
a linear dynamical system y = A p y that bounds the evolution of T, with A p defined as: 



A 






1 


. 


. 











1 . 


. 











. 


1 











. 


. 


1 


Xo 


Xi 


X2 ■ 


• Xp-2 


X P ~i 



(V.22) 



Using the trajectories of this linear system we rewrite the triggering condition Ji(t^, h p (z)) = 
as: 



yi(^VW) 



1 o 



A t-I- 

e p 



r(z) 



(V.23) 



It is clear that by fixing = in equation (V.23) we obtain an equation describing the set of 
states whose inter-execution times are lower bounded by t*. Hence the set flj t defined by: 

(V.24) 

represents a bound for the isochronous manifold Q tit in the following sense: 

t(z*) > t*. 

As mentioned before, we will use this bound to compute the inter-execution times through the 



Vli = {z* e R 2n : yi(t„/i p (^)) = 0} 



z * e fit 



scaling laws derived in section IV 



C. A self-triggering formula 



In order to apply the scaling law (IV6) with Cl^, we first find the intersection between the 
homogeneous rays and f^. Since we are searching for a self-triggering condition to be applied 
online, it is desirable to have a closed-form expression for those intersecting points. Towards 
this objective we state a simple lemma describing a useful property of the Lie derivative of a 
homogeneous function along homogeneous vector fields. 

Lemma V.3. Consider a map T : M — > M homogeneous of degree $ and a vector field 
Z : M — > TM homogeneous of degree £ G K. Then, the k-th Lie derivative of V along Z 
is homogeneous of degree $ + k^: 

C k z F(Xz) = \* +k!l £ k z r{z), z e M, A > 0. (V.25) 
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The result can be easily proven by induction. To make use of the scaling law (|IV.6[), we 



identify Xz with a point z* in and z with z(ti). To compute the intersection point between 
a ray passing through z(ij) and il^ we substitute the equation of the ray z* = Az(tj) in the 
triggering condition ( |V.23[ ): 

r(Az(t<)) 



yi(t*,/i p (z*)) = y 1 (t M /i p (Az(t i ))) 



1 0...0 



e Ap ** 



£ z r(Az(*i)) 
^-^(Az^)) 



p-i 



i=0 



with: 



Pi{U,z(U)) = e Ap S (i+ i) C* z T(z{U)), 



i = 0,...,p-l. 



(V.26) 



(V.27) 



Since we are interested in the positive real solution, we can divide equation (V.26) by A tf and 
define the variable q = A^ to obtain the equivalent equation: 



p-i 

i=0 



(V.28) 



The positive real solution of this equation can be substituted into the scaling law in (IV.6) to 
obtain a self-triggering condition, as summarized in the following theorem: 

Theorem V.4. Consider a dynamical system z = Z{z), homogeneous of degree £ > 0. Let 
U {oo} be the times implicitly defined by T(z(t(z), z)) = 0, where V is a homo- 



t : 



a 2n 



— >• 



geneous function. For any set of coefficients x ? Xii ■ ■ ■ > Xp-i £ ^ satisfying (V.12), and for any 
t* G M.q the submanifold fl^, defined by (V.24), upper bounds the isochronous manifold fit* in 
the following sense: 

z e t(z) > u. 

Moreover, the function : M — > M defined by: 
with q satisfying: 



(V.29) 



p-i 

Y,Pi(t*,z(ti))q l = 

i=0 



(V.30) 



provides a lower bound for the time r(z(i;)), with /^(i*, z(tj)) gz'ven fry (K27). 
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Equation ( |V.28[ ) reflects the tradeoff between complexity and accuracy of the computation 
of the inter-execution times: increasing p improves the approximation, but also increases the 
computational cost. The value p = 3 seems to be a sensible choice since it allows us to find a 
simple closed-form solution for ( |V.28| ). Equation < |V.29| ) then becomes: 

4_ 



(V.31) 



-ft + sign(/3 ) jft-lfafa 
Remark V.5. For non-homogeneous systems and non-homogeneous triggering conditions, a 



lower bound r^(z(tj)) for the inter-execution time r(z(tj)) is also given by (V.29), where now 
the coefficients ft are computed with the homogenized system and the homogenized triggering 
condition, as defined in section |7T] 

Remark V.6. In the case of linear systems, the inter- execution times remain constant along 
homogeneous rays since the degree of homogeneity is 0, as dictated by Theorem \IV.3 In order 



to obtain self-triggering formulas through Theorem V.4 the linear system has to be first rendered 



homogeneous of degree £ > 0, as explained in Section IV-B 



VI. Example 

We compare the results herein developed with our previous work in [[AT08J and HAT 1 OH . The 
equations of the example in [AT08| are: 



Xl 

x 2 



XiX 2 . 



Xl*2 



x 2 x 2 



(VI.l) 



with control law u 



X\x\. The operating region is a ball of radius 1 around the origin. 



We use the same triggering condition as in [AT08|, given by: 

|e| 2 = 0.0127V|x| 2 , <re]0,l[. (VI.2) 

Using SOSTOOLS [PPSP04], we obtain a set of coefficients x-i satisfying inequality ( |V13[ ) for 

p = 3: 

Xo = 105.970, xi = 0.021, X 2 = 1-033. 
These coefficients define the approximation of the isochronous manifold. Figure [T] depicts Qj f for 



£* = 1ms, computed according to equation (V.24), and the isochronous manifold Q tf computed 



via numerical simulations. As proven in Theorem V.4 the exact isochronous manifold encloses 
n^, since times enlarge as we approach the origin. Hence, two conclusions can be drawn from this 
figure. First, we notice that for the chosen value of t* the set nearly coincides with fi, while 
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a periodic 



self-triggered self-triggered 



lATlOl Thm. I V.4 1 



event-triggered 



0.1 0.39 0.39 1.50 1.55 

0.2 0.79 0.79 3.00 3.06 

0.3 1.18 1.18 4.50 4.58 



TABLE I 

Average time along a sphere for the BAT08B example (in ms.) 

there is a considerable gap between and the sphere, showing that the lower bound developed 
in HAT 10H was not tight. Moreover, the shape of the isochronous manifold is clearly different 
from the sphere, implying that even if a tight bound for the times on a sphere is computed, there 
exist many points where the inter-execution times will be significantly conservative. 

Since we are not aware of any other work addressing self-triggering strategies for nonlinear 



systems, in the simulations we compare the inter-execution times defined by Theorem V.4 against 
the periodic strategy, the event-triggered times generated by the triggering condition (III.7) and 
the self-triggering technique in HAT 1 OH . The period is computed as explained in [Tab07J. The 
trajectories are similar under the 4 different implementations for all tested initial conditions. 

As our previous work in HAT 1011 describes the exact evolution of the times along rays, we focus 
our comparison on the evolution of times across rays, which is the main topic addressed in this 
paper. For that purpose, we consider 20 initial conditions equally spaced along a sphere. The self- 
triggering condition developed in HAT 1011 defines the inter-execution times just as a function of 
the norm of the state, hence it will produce the same inter-execution times for all points lying on 



the sphere. On the other hand, Theorem V.4 takes into account more information contained in the 
state. Figure [2] shows the inter-execution times under the 4 strategies as a function of the position 
in the boundary of the sphere (that is, x = (cos(0),sin(0)), for 9 e [0, 2tt[). Table [I] represents 
the average inter-execution time for the 20 points along the boundary for different values of a 
{i.e., for different degrees of performance). We can observe both in Table [I] and Figure [2] how 
the new proposed technique nearly matches the event-triggered times and improves significantly 
the times generated by the former self-triggering strategy, which generates in this case the same 
times as the periodic strategy. 

Finally, Figure [3] shows the evolution of the inter-execution times along trajectories for the 4 
different implementations for a particular initial condition (x = (0.4, 0.7)) and for a simulation 
of 5s. Again, we can observe that the new self-triggering technique nearly tracks the evolution 
of the event-triggered times. 
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Fig. 1. Reference manifold for IAT10I . isochronous manifold and its approximation according to equation 
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Fig. 2. Inter-execution times along a sphere of unitary radius. 

VII. Accuracy of the self-triggering technique 

The derived self-triggering formulas represent an accurate bound for the event-triggered times 
in the example in Section |VI| However, there is no guarantee that such bounds will be accurate in 



general. Lemma V.2 provides a constructive procedure to derive a linear model lower bounding 
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time(s.) 
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Fig. 3. Evolution of the inter-execution times along trajectories. 

the evolution of the triggering condition. The same procedure can be used to derive an upper 
bound t^, by reversing the inequality ( |V. 12[ ): 

p-i 

(C p z T)(z,w) >J2xi(^)(z,w). (VII.l) 

i=0 

Using this upper bound r T , a bound for the mismatch between the self-triggered times and the 
event-triggered times can thus be computed: 



< r(z(t,)) - t^z(U)) < r^zfe)) - t\z(U)). 



(VII.2) 



While the self-triggered times generated by Theorem V.4 for p < 3 might be useful for some 
applications, they could be conservative in other cases. As mentioned before, self-triggered 
formulas are expected to yield larger inter-execution times as the value of p is increased, that 



is, when high order approximations are used in (V.13 ). However, for p > 3 no simple closed 



form expression can be obtained for equation ( |V.28[ ). To overcome this drawback, in the next 
subsection we develop an iterative algorithm that considers high order approximations (p > 3) 
to compute a self-triggered formula. 



A. An iterative algorithm 



To compute the next inter-execution time r(z(tj)) equation (V.28) needs to be solved online 



upon the measurement of the augmented state z(tj). Equation ( |V.28[ ) is polynomial of order 
p — 1 in the variable q. For p > 4, closed-form expressions for r are too complicated to be 
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evaluated online. There is a vast literature on algorithms that find approximations of the roots of 
polynomials. However, we seek approximate solutions that are always lower bounding the exact 



root of (V.28), in order to obtain inter-execution times that guarantee stability of the system. 
For that matter, we develop an iterative algorithm that always produces lower bounds for the 
inter-execution times. 

For convenience we define the map: 



S x (x) := (xi, Ax 2 , 



X p x r 



(VII.3) 



Notice that Lemma V.2 allows us to derive bounding models of different orders. We first 



construct two bounding models using Lemma V.2 a low order model y, with order p < 4, and 
a high order model y, with order p > 4: 



e A ^/I(z) 



(VII.4) 



The low order model will be used to compute the solutions and the high order model will serve 
to evaluate the accuracy of the solution, in the spirit of the Newton-Raphson method. 



T(z(t*)) 




T(t,z(ti)) 



Fig. 4. Iterative algorithm for the computation of the inter-execution times. 

For clarity of exposition, we first describe the main ideas behind the algorithm. The procedure 
is sketched in Figure |4j Both y i and y 1 upper bound the evolution of the exact evolution of the 
triggering condition T, hence their roots lower bound the exact solution r(z(ti)). The first bound 
for the inter-execution time can be computed directly using the low order model: 

rl = minarg^y^/i^))) = 0) (VII.5) 
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This value represents a bound on the inter-execution time r(z(tj)) defined by the triggering 
condition, r(z(r(z(ij)), z(U))) = 0. Since y does not represent in general a tight bound for the 
evolution of T, it is expected that the triggering condition has not been satisfied at rf. To refine 



the self-triggered time in ( |VII.5[ ), we evaluate the high order model at this instant of time: 

«i = y(7^(zfe))) (VII.6) 
The high order model provides us with a lower bound for the gap between the exact inter- 



execution time as defined by the triggering condition and the first value obtained in QVII.5| ) 
through the low order model. We use this value v\ as the new initial condition for y to improve 
the initial guess rj: 

r*(zfe)) = rf + rt >rj = minarg^t, Cv x ) = 0) (VII.7) 

with C = [ Ii xp 0i x (p_p) ]. The procedure can be repeated iteratively to obtain more accurate 
self-triggered times, as stated in the following theorem. 

Theorem VII.l. Consider a dynamical system z = Z{z), homogeneous of degree £ > 0. The 
time r^(z(ti)) given by: 

T^z(t l )) = f^rt (VII.8) 



i=i 



(VII.9) 



with rf = minarg r (yJr,Cvi- X ) = o) , C= J lx? lx( p_p) 
and v = fi p (z(U)) (VII.10) 
v r = y(7^,« P _i), r = l,2,...,n der -l (VII.l 1) 

represents for any n iter > a lower bound for the inter- execution times r : M. 2n — > IR + U {oo} 
implicitly given by the homogeneous triggering condition r(z(r(z(tj)), z(tj))) = 0, with y and 
y as defined in ( |VII.4| ). 



Proof: According to Lemma V.2, y represents a bound for the evolution of the triggering 
condition: 

r(z(*,z(*i))) ^y&vtyzfa))). (vn.i2) 

Hence we can use y to compute a first bound for the inter-execution time r(z(tj)): 

ri := minarg^y^t,^^^))) = 0) < r(z(t i )) := min arg t (r(z(t, z^))) = 0). (VH.13) 
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This bound will be refined using the high order model. According to Lemma |V2| , y also represents 

(VII. 14) 



a bound for the evolution of the triggering condition and its Lie derivatives: 



C l ^T(z(t, zfo))) < y,(*, ti p {z(U))), 1 = 1,..., p. 



We define 



^:=y(Tj",//(z(i,))) (VII.15) 
which represents a bound for F and its first p derivatives at the time instant rf. We shall use 



this bound as the new initial condition for y . We now rewrite inequality ( |VII. 13[ ) using v\\ 

4 : = minarg^y^, CvJ = 0) < r 2 := min arg t (r(z(t, Zl )) = 0) (VII. 16) 

where z\ represents the extended state after flowing along the vector field Z for rf units of time: 

Zl = z(rf,z(U)). (VII.17) 



Notice that inequality ( |VIL16[ ) holds since Lemma V.2 also applies whenever the initial conditions 
satisfy: 

w-i.v. , ^ .. - (VII. 18) 



C t z 1 T(z )<y 0l , 1 = 1,..., p. 



which is the case in (VII. 16) since: 



C l ^T( Zl )<v h , 



(VII. 19) 



Merging equations ( |VII. 13[ ) and ( |VII. 16| ) we obtain a second bound for the inter-execution time: 

rf + 4 < rf + r 2 = t(z(U)). (VII.20) 

Applying the procedure iteratively we finally obtain the final expression in ( |VIL8[ ), with rf = 
minarg^y^tjCvi-i) = 0). ■ 

Remark VII.2. There is no guarantee that y < y, for there does not exist a partial order on the 
approximations (that is, a high order approximation does not need to be more accurate than a 
low order approximation for all points in the extended state space). Hence, the bound does not 
need to increase monotonically with respect to the number of iterations. However, the iterative 
algorithm can be easily modified to account for such possibility: if > at some step i, the 
computation should stop at this point, as the high order model estimates that the triggering 
condition has already been satisfied and thus will yield shorter times. 

In practice, since the expression for y is transcendental in t, a closed-form expression for rf 
and T2 cannot be directly computed, as explained in Section V-B To find closed-form expressions 
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for such times, we proceed as in Section fV-C] we use the scaling properties of the inter-execution 
times for a homogeneous system. The approximation of the isochronous manifold f2 using y is 
defined by: 

y i (U,^(z*))=0. (VH.21) 
We find the intersection of the isochronous manifold with the homogeneous ray = Az(tj): 

y^^XziU))) = yfi^S^Coo)) = 0, o = (jF(z(U)) (VH.22) 



since fjP-(\z(ti)) = 5 x (fi?-(z(ti))), according to Lemma V.3 Equation ( |VII.22[ ) is now polynomial 



of order p in A (after manipulating terms as in Section V-C ), and therefore a closed-form solution 
exists (since p < 3). We denote such solution as Ai: 

Ax = arg^y^c^Coo)) = 0) < arg A (z(*„ Az(t t )) = 0). (VH.23) 

Equation ( |VIL23| ) implies that the point Aiz(tj) has an inter-execution time bounded by t*: 

r(A 1 z(t i )) > (VII.24) 



And using the scaling laws from Theorem IV.3 we compute the first bound for the inter-execution 
time at z(ij): 

r(z(ti)) >t 1 = \{t* (VII.25) 



which corresponds to equation ( |V.29[ ), with q = X\. Hence, for n. 



iter 



1 we obtain the bound 



in Theorem V.4 Likewise, we can use the same technique to compute the second iteration. As 
before, we define: 



oi =y(t*,5 Xl (o )) 



(VII.26) 



which represents a bound for T and its p Lie derivatives at the time instant t*. We rewrite the 

(VII.27) 



left hand side of ( |VII.22[ ) using o\\ 



Tl {U,5\C 0l )) 



0. 



We solve for A as before: 



A 2 = arg A (y l (t, ) 5 A (Co 1 )) = 0) < arg A (z(t„ Xz a ) = 0) 



(VH.28) 



where z a represents the extended state after flowing along the vector field Z for £* units of time, 
starting with initial condition Aiz(tj): 



z a = z(U, Aiz(ij)). 



(VII.29) 
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Equation ( |VII28| ) implies that the point X 2 z a has an inter-execution time bounded by t*: 

r{X 2 z a ) > U. (VII.30) 

Resorting to the scaling properties for r again we can compute a bound for the inter-execution 
time at z(tj): 

T(X 2 z a ) = Afr(z ) = A ? 2 r(z(t, X&fc))) = A 5 2 (r(A 1 z(t i )) - *♦) = A|(A^r(z(t i )) - Q > t*. 

(vn.3i) 

Finally, solving for r(z(tj)) we obtain: 

t(z(U) > X{U + (VH.32) 



Applying the procedure iteratively we obtain the following corollary, analogous to Theorem VII. 1 



Corollary VII.3. Consider a dynamical system z = Z{z), homogeneous of degree £ > 0. The 
time T^(z(ti)) given by: 

( niter / I \ \ 

e (n M j j *. ( vn - 33 ) 



A j = minarg A (y i (t„5 A (Co J _ 1 ))=0), (VII.34) 
an J o = fJ?{z{ti)) (VII.35) 
o r = y(U,6 Xr (o r ^)), r = 1,2,..., 71^.-1. (VII.36) 

represents for any m ter > a lower bound for the inter-execution time r : IR 2n — > IR + U {oo} 
implicitly given by the homogeneous triggering condition r(z(r(z(tj)), z = 0, wzY/z y and 
|7 a* defined in ( VII.4| ). 



Remark VII.4. 77ze proposed algorithm has the property of being anytime, since it can be 
stopped at any iteration and still provide a bound for the inter-execution times. Other numerical 
algorithms, such as any Householder's method l \Hou70V . could be used to compute an approx- 



imation X for the solution of the polynomial equation (V28); then, the inequality Vi(X) < 
needs to be checked to guarantee that the solution indeed bounds the exact solution. However, 
in this approach no conclusions can be reached should the inequality not be satisfied. 
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B. Example 

We consider now the control of a rigid body, previously studied in HAT 1 PL The state space 
representation of such system with two inputs can be simplified to the form: 



Xl 

x 2 

X.3 



Ui 

u 2 

XlX 2 . 



(VII.37) 



The following nonlinear feedback law was designed in HBI89M to render the system globally 
asymptotically stable: 



— x\x 2 - 2x 2 x 3 — xx — x 3 

„2 



Ui 

u 2 = 2xiX 2 x 3 + 3x1 ~ x "i- (VH.38) 
We use the same triggering condition as in HAT 1 PL given by: 

(VII.39) 

Since the system is not homogeneous, we introduce an auxiliary variable w, as explained in 



0.79V|x| 2 , <re]0,l[. 



Section IV-B The homogenized closed loop system becomes: 



Z(z) 



-(xi + e x )(x 2 + e 2 )w - 2(x 2 + e 2 )(x 3 + e 3 )w - (x x + 



w 



(x 3 + e 3 )w 2 



2(xx + ei)(x 2 + e 2 )(x 3 + e 3 ) + 3(x 3 + e 3 ) 2 w - (x 2 + e 2 )w 2 

XiX 2 W 

(xi + ei)(x 2 + e 2 )w + 2(x 2 + e 2 )(x 3 + e 3 )w + (xi + ei)w 2 + (x 3 + e 3 )w 2 
-2(xi + ei)(x 2 + e 2 )(x 3 + e 3 ) - 3(x 3 + e 3 ) 2 w + (x 2 + e 2 )w 2 

— XiX 2 W 



with z = (x, e, w) representing the extended state. It is not necessary to homogenize the triggering 
condition since it is already homogeneous. 

We analyze the accuracy of the self-trigger times with respect to the number of iterations 
executed. The degrees of the low order and high order models are p = 3 and p = 4 respectively. 



Using this augmented vector field Z and the triggering condition ( |VIL39| ) we obtain a set of 
coefficients satisfying ( |V.13| ) for p = p and p = p: 



P = P- Xo = -73.2528, Xl = 1.7157, xa = 1-8299. 
p = p : Xo = -57.8151, xi = 1-4923, X2 = 24.9920, Xs = 47.6313. 
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niter = 1 


self-triggered 

T^iter — 2 


niter = 3 


pvpnt-tn CT(*prpfl 

L v L 1 1 I 11 1 L,!- 1 - 1 L LI 


0.50 


216.08 


220.14 


220.18 


220.34 


0.65 


267.53 


278.29 


278.81 


285.41 


0.80 


322.12 


349.97 


354.88 


355.75 







TABLE II 

Average time for the rigid body example (in ms.). 



With this set of coefficients we can finally construct the polynomial equation in ( |V.30| ), with 
zfc) = (x(^),e(£j), 1) (since w(t) = 1). 

The proposed self-triggering technique is evaluated for 25 initial conditions equally spaced 
along the boundary of a sphere of radius 1. Table [XT] shows the inter-execution times (averaged 
over the considered initial conditions) when the number of iterations ranges from 1 to 3. It 
can be seen how the mismatch between the self-triggered times and the event-triggered times is 
quickly reduced as the number of iterations grows. 

VIII. Selection of design parameters 

To conclude the paper, in this section we study the role played by the different design 
parameters of the proposed procedure. 



A. Choice of the embedding degree 



In Section IV an embedding procedure was derived to render a system and a triggering con- 
dition homogeneous of any desired degree. In practice, one would like to choose the appropriate 
degree of the embedding to obtain accurate bounds for the inter-execution times. It was shown in 



Section IV-C that the scaling law does not depend on the degree of homogeneity of the triggering 
condition. In the following theorem we prove that embeddings of different degrees generate the 
same bounds for the inter-execution times of the original system. 



Theorem VIII. 1. Let n : 



i>2n+l 
Z 

w 



— > 



U {oo} be the inter-execution times of the system: 



Si(z,w) 



w* 1+1 0(w -1 z) 



(VIII. 1) 



implicitly defined by the homogeneous triggering condition T(w 1 z(r 1 (^ , w ),z , w )) = 0, and 
let r 2 : 



»2n+l 



U {oo} be the inter- execution times for the system: 

76- 



z 
w 



5*2(2, w) 







(VIII.2) 
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implicitly defined by r(w -1 z(r 2 (.2o> wq), z , w )) = 0. Consider two bounds rf and r 2 for t\ 
and r 2 , computed according to Theorem V.4 and valid in the sets Q\ and Q 2 = {(z,w) G 
]R 2ri+1 |(u7^ 2 ^ 1_1 z, w^ 2 ^ 1 ) G Qi} respectively. Then, the two bounds coincide for w = I: 

rf( Zo , 1) = 4(z , 1), v(2 , 1) e Oi n n 2 . (vni.3) 

Proof: Since the degree of homogeneity of the triggering condition T does not play any 



role in the inter-execution times (as explained in Section |TV-C[ ), we assume it to be zero, without 
loss of generality. The bounds r| and r 2 are given by: 

QjU, J = 1,2 (VIII.4) 



ECo = 0, frfa, 1) = e^\ (i+1) (C s T)(z , 1), % = 0, . . . ,p - 1. 



A P2 and the Lie derivatives of 



Thus both bounds will coincide if fi ix = fy 2 , that is, if A Pl 
both vector fields coincide for w — 1, as imposed by the following equations: 

(c Sl r)(z , 1) = (c S2 r)(z , 1) i = 0, 1, . . . P - 1 

(£^r)(^ 1)Wl ) < ^^(^H^i^i) V^i.^O g fix 

i=0 

p-1 

^(£5 a r)(z2,w 2 ) < X)xi(4 a r)(22,w 2 ) v(z 2 ,w 2 ) g fi 2 

i=0 

where the coefficients Xi are the same for both systems (in order to have the same A p matrix). 



(VIII.5) 



(VIII.6) 



To prove equality ( |VIII.5[ ), we notice that the submanifold defined by w = 1 is invariant 
under the flows of Si and S 2 . Since Si and S 2 coincide for w = 1, the evolution of the 
triggering condition in this submanifold is the same under both vector fields, and therefore the 
Lie derivatives coincide for any (z, 1). 

To prove the double implication ( |VIII.6[ ), we resort to the homogeneity properties of the Lie 



derivatives of T, as stated in Lemma V.3 The proof relies on the fact that the vector fields 
coincide for w = 1 and at the same time they are determined by their restriction to w = 1. 
To proceed, we represent a point (zi,Wi) G Oi as Xi(z a , 1), and we rewrite the left hand side 
of ( [VIII.6D as: 



p-i 



jC p s T(Xiz a ,\i) < ^x*£kr(Ai;z ,Ai) VAi > 0|(A^ o ,Ai) G Qi. 



(VIII.7) 



i=0 



Since the z-th Lie derivative £* Sl T is homogeneous of degree we can rewrite ( |VIIL7[ ) as: 

p-i 

X p ^C p s r{z a ,l) < ^Xi^ 1 4^(^,1) VAi > 0|(Ai* a ,Ai) G fii. (VIII.8) 



i=0 
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Defining A 2 := A^ 1 ^ 2 and using equality ( |VIII.5[ ) we obtain: 

p-i 

X p ^C p s T(z a , 1)<J2 ^fC S2 T(z a , 1) VA 2 > 0\(Xi 2 ^z a , Aj /€l ) G Q 1 . (VIII.9) 



i=0 



Using the homogeneity properties again for the Lie derivatives equation (VIII.9) becomes: 



jC p s T(X 2 z a ,X 2 ) < ^xAr(k,A 2 ) VA 2 > 0|(A? /a z a ,Af /6 ) G fix. (VIII.10) 



i=0 



The point (z 2 ,w 2 ) := (X 2 z a , A 2 ) is related to {z\,Wi) according to: 



(z 2 , w 2 ) = A 2 (z a , 1) = Xf 



Ai 



zi, 1 



A x ^1,A 1 



wf^zuwf 1 ^ 



(VIII.l 1) 



Thus, inequality ( |VIIL10D holds for all («i,twi) G fli, that is, for all (wf^ 1 K - 



z 2 , u> 2 "' ") G Qi. 

Hence, we can conclude that the left hand side of ( |VIII.6| ) for (zi, Wi) G f2i implies the right hand 



side of ( |VIIL6[ ) for (-2 2 ,w; 2 ) G Vt 2 . The other side of the implication can be proved analogously 
by swapping the subindices 1 and 2. Therefore, r| and coincide for all those points (z , 1) 
where the bounds are valid, i.e., for all (z , 1) G fix fl £l 2 ■ 
The set f2 2 can be thought as an inflation of the set f2i as dictated by the different degrees of 
homogeneity. As both r| and coincide for all the points satisfying w — 1, a bound for the inter- 
execution times of the dynamical system z = Z{z) under the triggering condition T(z(r, zq)) = 
can be equally computed using either system ( |VIIL1[ ) or ( |VIII.2[ ). In other words, the degree of 
the homogenized system does not affect the accuracy of the bound for the inter-execution times. 



B. Choice oft* 

In order to derive the self-triggering technique, we used the concept of isochronous manifolds 
in Section |V} It was proven in Proposition |V1 that such manifolds exist for any value of £*. 
This parameter has to be selected offline since computing e Apt * online would imply a high 
computational cost. Should the exact computation of this manifold be possible, any value of 

would yield the exact times. However, such exact computation is in general not feasible. 



Lemma V.2 describes a way to compute one approximation bounding the exact isochronous 
manifold. Hence, a natural question is whether the choice of t* plays a role in the accuracy of 
the times, and, if so, which value of generates the largest times (and thus the most accurate 
bound). 



To proceed, we analyze the role played by in equation ( |V28[ ): 

p—i p— i 

Y.&MUM = 5> A ^ 1(m) C z T{z{U))q* = 0. 



(VIII. 12) 



i=0 



i=0 
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Taking into consideration the special structure of the A p matrix (as defined in ( |V.22[ )), we can 
write the Taylor expressions for the elements of e Apt * that appear in ( |VIII. 12[ ): 



oo 

e ApU i(i + D = E Tfi( A iUi + D = -fi + 0{t 



(VIII. 13) 



where 0(t%) denotes terms that depend on n-powers of t*, for n > p. Using ( |VIII.13| ), equa- 
tion ( |VIIL12[ ) becomes: 

£ rfl C z T(z(t t ))q* + Y^O(tOC z T(z(U))q^ = 0. (VIII.14) 

i=0 ' i=0 

Since we are interested in the dependence of on i*, we rewrite equation (|VIII.14|) using 



qt* 



p-i 



E l^y &n<u)) + E°(^)Ar(z(t,)) 

i=0 ' i=0 

E k^y c z v{*{u)) + ^0(04^))^ 



o. 



(VIII. 15) 



i=0 i=0 

The first term of the previous summation corresponds to the first p terms of the Taylor series of 
the triggering condition T. Therefore, the self-triggering formula that has been proposed in this 
paper represents a modified Taylor series that guarantees a lower bound for the inter-execution 
times. Likewise, we can obtain a similar expression for an upper bound r^, as mentioned at the 



beginning of Section VII 



E c z T Wi)) + P ^o{tr)c z T{z{t l )){^y = o. (vm.i6) 



i=0 



i=0 



Since the first terms of ( |VIII. 15[ ) and ( |VIII. 16[ ) coincide, the gap between the expressions for the 
upper bound and the lower bound is a function of 0(t*), and therefore the mismatch between 
and -H- increases with t*. However, cannot be chosen as small as desired. The set flf^ 



defined in (jVTTj) represents a bound for the isochronous manifold, valid in the set M where 
the inequality ( |V13[ ) holds. As decreases, the diameter of the set flj t increases, since times 



enlarge as we move far from the origin (as dictated by the scaling law ( |IV6[ )). Hence, has to 
be chosen so that is contained in the operating region M. Therefore, a sensible choice for 
U is as follows: 



t* = minjt G 



flj C M}. 



(VIII. 17) 



Since the computation of ( |VIIL17[ ) can be difficult in general, in practice t* can be chosen as the 
minimum value for which flj C M can be guaranteed. Notice that the choice of = is not 
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a possibility, since there does not exist a set of coefficients \i sucn tnat (|V-13|) holds globally. 



Indeed, the left hand side of the inequality ( |V.13[ ) is a homogeneous function of higher degree 



than the homogeneous functions on the right hand side, thus there always exist points where 
C P Z T is greater than any linear combination of the previous Lie derivatives. 

IX. Conclusions 

In this paper we developed a new self-triggering technique for nonlinear systems. Our previous 
work was extended in two directions: first, we enlarged the class of systems and triggering 
conditions that can be considered. Second, the new self-triggered formula reduces the conser- 
vativeness of previous techniques by resorting to the concept of isochronous manifolds, herein 
introduced. These results help consolidating self-triggered control as an interesting alternative 
to event-triggered control whenever the latter is not implementable. Moreover, even in those 
implementations where event-triggered control is preferred (e.g. for robustness purposes), the 
self-triggering approach provides the tools to analyze the time requirements imposed by event- 
triggered control. 
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